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Abstract 



The Raman interaction of a trapped ultracold ion with two travelling 
wave lasers is studied analytically with series solutions, in the absence of the 
rotating wave approximation (RWA) and the restriction of both the Lamb- 
Dicke limit and the weak excitation regime. The comparison is made between 
our solutions and those under the RWA to demonstrate the validity region of 
the RWA. As a practical example, the preparation of Schrodinger-cat states 
with our solutions is proposed beyond the weak excitation regime. 
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I. INTRODUCTION 



The preparation of ultracold ions as well as the production of nonclassical motional states 
of these ions plays a central role in ion trap experiments^ 1 ^ . As the trapped ion system is 
in nearly perfect isolation, immune from the violation of the environment, and meanwhile 
the electric-dipole forbidden transition is usually introduced for avoiding the dissipation of 
the ion, the decoherence in the preparation and preservation of nonclassical states can be 
effectively suppressed. Consequently, the ion trap system is also being considered to be a 
promising tool for the quantum computing' 4-6 ' . 

The various theoretical schemes for generating the nonclassical states of motion of the 
trapped ions and achieving quantum computing with trapped ions are based on the Jaynes- 
Cummings (JC) model' 7 !, in which the ion is simplified to be of two levels, the trap's potential 
is quantized as a harmonic oscillator, and the radiating lasers are supposed to be classical 
forms of standing or travelling waves. The general consideration is taken for the case of 
Lamb-Dicke limit (LDL) under the weak excitation regime (WER), which corresponds to 
the actual case in present ion trap experiments' 3 ' . The LDL means the ions moving within 
the region of the space much smaller than the effective wavelength of the laser wave, and 
the WER means the laser-ion interaction much smaller than the trap frequency. For this 
case, some techniques developed in the framework of cavity QED based on JC model can 
be immediately transcribed to the ion trap system by taking advantage of the analogy be- 
tween the cavity QED and ion trap problem. In the case of motion of ions exceeding the 
LDL, we can use the technique proposed in Ref.[8] for a treatment. In a word, under the 
rotating- wave approximation (RWA), the problem of trapped ions interacting with lasers is 
solvable. However, the strength of the coupling between the ions and lasers can be conve- 
niently adjusted simply by changing the intensity of lasers. If the coupling strength is much 
larger than the trap frequency, called the strong excitation regime (SER), then the RWA 
is no longer applicable since the rapidly oscillating (i.e., counter- rotating) terms also make 
significant contribution in the interaction, which is much different from the case in Cavity 
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QEDl 9 !. It has been reported that the SER is useful for preparing nonclassical motional 
states of a trapped ion and implementing quantum computing with trapped ions! 6 ' 10-12 '. 
For example in Ref.[10], Schrodinger cat states could be prepared readily in simple way 
in the SER. Based on that work, a proposal was made for motional state engineering and 
endoscopy in the SER^ 11 '. As the quantum state is sensitive to decoherence, the manipula- 
tion of trapped ions in the SER is of importance to the experimental implementation and 
measurement because of the great reduction of the operation time. Due to the same reason, 
a scheme for the hot-ion quantum computing also requires the condition of the SER' 6 !. In 
contrast to the theoretical studies, the ion-laser interaction in the SER has not yet been 
reported experimentally since the high-laser intensity produces increased off-resonant exci- 
tation of carrier transitions which limits the final ground state occupation of the cooling 
process' 13 ^, and the strong laser also modifies the original electronic levels of the ions' 14 ', 
making the problem somewhat complicated. In a numerical study of the strong Raman 
sideband excitation driving the vibrational ground state of the trapped ion' 15 !, it is shown 
that the non-RWA interaction terms produce some limitations for two strong laser beams 
detuned by the vibrational frequency to be used as Raman motional displacement beams, 
whereas the interaction is useful for the preparation of the Schrodinger cat state. 

Generally speaking, as the RWA is invalid, the ion trap problem in the SER can not 
be solved analytically with JC model. So some approximations have to be introduced in 
Refs. [10,11] and the studies in Refs. [12,15] are resorted to numerical calculations. In fact, 
as shown in Ref.[16], even if the ions are governed in the WER, in which the RWA is valid, 
the deletion of the RWA can also present some interesting results much different from those 
under the RWA. In that work, some particular analytical solutions have been obtained in 
the absence of the RWA for large detuning cases, and Schrodinger cat states could be readily 
prepared. However, that work is too simple. We hope to investigate the Raman configuration 
which has been extensively applied in actual ion trap experiments' 3 ', instead of the simple 
case in Ref.[16]. We also hope to exclude the assumption related to both the WER and 
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LDL. Therefore, in this contribution, we will try to undertake a more general and complete 
investigation on the trapped ion under the Raman process. Moreover, a comparison of our 
results with those under the RWA will be made. Besides, we will propose a scheme for 
preparing the Schrodinger cat state beyond the WER. 

II. MODEL AND SOLUTION 

We study the single ultracold ion radiated by lasers in Raman-A-type configuration, 
as shown in Fig.l. The electronic structure is employed with two lower levels |e > and 
\g > coupled to a common upper state \r >, and the two lasers with frequencies uj\ and 
uji respectively are assumed to propagate along opposite directions. For a sufficiently large 
detuning to the level \r >, \r > may be adiabatically eliminated, and what we have to treat is 
an effective two-level system, in which the lasers drive the electric-dipole forbidden transition 
| g >^ \e >. The dimensionless Hamiltonian of such a system in the frame rotating with 
the effective laser frequency coi(= uj\ — UJ2) can be written as' 10 ! 

H =^-a z + a + a + ^(a + e ir <* + a_e~^) (1) 

where the detuning A = (u — iO\)jv with u being the optical resonance frequency, i.e., 
the transition frequency of two levels of the ion, and v the frequency of the trap. Q is 
the dimensionless Rabi frequency and r\ the effective Lamb-Dicke parameter given by r\ = 
rji + r] 2 with subscripts denoting the counterpropagating laser fields. <7j (i = ±, z) are Pauli 
operators, x = a + + a is the dimensionless position operator of the ion with a + and a being 
operators of creation and annihilation of the phonon field, respectively. The notations '+' 
and '-' in front of irjx indicate the absorption of a photon from one beam followed by emission 
into the other beam and vice versa, respectively, v is generally supposed to be much greater 
than the atomic decay rate, called the strong confinement limit, for neglecting the effect 
of the atomic decay. In general, one may expand e ±irix in Eq.(l) to the first-order terms 
of rjx and neglect other higher-order terms by supposing the ion governed under the WER 
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(Q <C 1) and within the LDL (rj <C 1). However, as we hope to treat the problem more 
exactly and generally, we first perform following unitary transformations on Eq.(l), that is 

H 1 = V + U + HUV = -a z + a + a + g(a + + a)(<7+ + <7_) + e(<r+ + cr_) + g 2 (2) 
( D{(5) -D(/3)\ 

where U = \ and V = -4- e -*™ +a / 2 with £>(/3) = e *K« + +<0/2 g = r}/2 

and e = —A/2. Actually, the unitary transformation [/ made above is identical to that 
in [17]. But for coinciding with the standard form of non-RWA model, V is performed. 
Comparing with the standard non-RWA JC model' 9 ', there exist an additional driven term 
and an additional constant term in Eq.(2), and the optical resonance frequency is replaced by 
the Rabi frequency. In what follows, we will adopt the coherent state representation' 18 ] , by 
which the non-RWA JC models have been analytically treated' 16,19 ]. Thus Eq.(2) is rewritten 
as 

H = % z + a^- + g(a + -^-){<J+ + <7_) + e(a+ + cr_) + g 2 (3) 
z da da 

where a is a complex number, and / da ^f exp(— |o:| 2 )|o:* >< a*\ = 1 with the relation 
between the coherent state and Fock state 

r dada* , , ,,,, » 1 _ 
n>= I — exp(-\a\ 2 )\a* > ■ ( 4 ) 
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Using the same idea as in Refs. [16,19], we assume the eigenfunction of the Schrodinger 
equation of Eq.(3) to be the series form 

/Vi(a) \ ( exp(-za)Y,n=o h nU 

V^ 2 (a); V e ^(-^a)En=oCna n , 
where z, b n and c n are constants determined later. Taking ty(a) into the Schrodinger equa- 
tion of Eq. (3) will yield two recurrence relations 

b n +i = ^— ^ ^ [(E + ^ - n - g 2 )c n + (gz - e)b n - gb n -i + zc_i], (5) 

Cn+i = 7 1 tt [(-E - 77 - n ~ 9 2 )K + (gz - e)c n - gc n _ x + z6 n _i] (6) 
g{n + l) 2 



where E is the trial solution of the eigenenergy of Eq.(3). As the procedure of solution is 
similar to that in Ref.[16], in what follows, we just list the results: 

(1) b n = c n = for n > 2, we have two solutions. One is E\ = 1 + e with E\ the 
eigenenergy of the system in this case, corresponding to z — g and a restricted relation 
of f2 = 2y/l + 2e — 4g' 2 . The other is E[ — 1 — e, with z = —g and a restricted relation 

of Q = 2y/\ — 2e — 4g 2 . The eigenfunctions for above two cases are respectively ty(a) = 

(bo + b±a\ (b + b 1 a\ 

and e flQ , , , where c = \^Xt-%^ bl = % l + f ~ 

+ e)c + (£ 2 - e)6o], c = ~ ^Z^X and &; = ±[(1 - e - + f )c - (^ 2 + e)6 ]. b 
and 6q can be set to be 1 for normalisation. 

(2) b n = c n = for n > 3, the two solutions are E 2 = 2 ± e corresponding to z — ±g 
and two independent restricted conditions, which can be solved from the following equation, 
respectively 

2g 2 - (1 - 2g 2 ± 2e)(2 + Sl/2 - 2g 2 ± 2e) + SI + fi 2 /4 
2c/ 2 - (1 - 2# 2 ± 2e)(2 - Q/2 - 2g 2 ± 2e) - ft + fi 2 /4 ~ 

V(2 + 3fi/4 - 2g 2 ± 2e) - VL(g 2 =f e)(3 ± 2e + ft - 2# 2 ) 



4# 2 (2 - 3fi/4 - 2# 2 ± 2e) + H[(l + H/2)(2 - ft/2) + 2(# 2 q= e) 2 - 3(# 2 =f e)] 

The eigenfunctions at this moment in the unit of bo = b' = 1 are, respectively, 

1 + b x a + b 2 a 2 



(7) 



(a) = e^ a 

and 



Co + CiGC + &2Q; 2 

1 + c^ck + 6 2 a; 2 



c + c x ot — b 2 a 2 

where cn - ( 1 -V+^)( 2 -^/ 2 - 2 g 2 +^)+^"" 2 / 4 -V Cl - lr 2 + e _ 2 _ n , / 2 _ ^ 1 6 _ 
wnere c - ( i_2 9 2 +2e )(2+fV2-29 2 +2e)-n-n 2 /4-2 9 2 ' c i — 9 L z + e 9 2 + v# ej^oj, 01 - 

i[(2 + e -^ 2 + f)co + ^ 2 - e ],62 = 2^{(^-e)(3-fi + 2 e -2^ 2 )+^ 2 + [(l + e -^ 2 ) 2 -l-f - 

^-n 2 -^ 2 -^^ 2 ^ \ r — 29 2 -CL-2g 2 -2e)(2-n/2-2g 2 -2e)-n+n 2 /4 ' _ 1 [o _ „2 _ £2 _ / 2 , \ ' 1 
2 9 \9 t ) J c 0i, C — ( 1 _2 ff 2_2 £ )(2+fi/2-2 a 2 -2e)-fi-fi 2 /4-2s 2 ' c l ~ gl^ t 9 2 \9^ t ) c Q\-, 

6' 1 = i[(2- e -^ 2 + f)c -^ 2 + e ],and b' 2 = ^{(g 2 + e)(3 - ft - 2e - 2^ 2 ) + g 2 + [(1 - e - 

g 2 ) 2_ 1 _V_n_ g 2_ {g 2 + e) 2 ]c > o} . 

(3) when above technique is extended to the case of b n — c n — for n > N + 1, the 
series solution becomes very complicated. However the eigenenergy of the system is still 
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of the simple form, that is, En — N ± e, corresponding to z — ±g and two independent 
sets of much complicated expressions of restricted conditions which we have to omit here, 
respectively. 



III. DISCUSSION 



As limitations for both the Rabi frequency and Lamb-Dicke parameter as well as the 
RWA are excluded in the present treatment, the results we obtained above are more exact 
and general than those depending on the RWA or the approximate expansion of small 
values of i] and Q. However, as the problem involving the counter-rotating terms is non- 
integrable' 16 ' 18 ' 19 ', the present solutions obtained by the termination of the series of the trial 
solution ^(a) are some particular ones, along with some restricted conditions related to 
some parameters of the system. Nevertheless, with the present solutions, we can investigate 
some cases in a wide range of parameters, particularly for the case beyond the WER and 
LDL. It is also of interest to make a comparison of our results with those under the RWA. 
Let us first consider a special case with both ~ e > 1 and rj — > 0. In such a case, Eq.(3) 
is reduced to 



Repeating the above procedure of the series solution for the simplest case, i.e., b n = c n = 



with the solution in a Fock state representation for n = 1. However, to make a comparison 
for eigenfunctions with former solutions, we have to return to the original representation 
before Eq.(2). The eigenfunction in the original representation is 



H s = -a z + a— + e(<7+ + <r_). 




and E s = 1 ± ye 2 + We can easily find that the eigenenergy E s is in good agreement 
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where we have used the fact' 16 ' 18 ] that e ga a n in a coherent state representation corresponds 
to the displaced Fock states \n, g > by means of Eq.(4) and Baker-Campbell-Hausdorff 
formula. Therefore, in the case of Q <C 1 but e ^> 1, we have from Eq.(8) ^ os ~ or 

w 

, which means no transitions existing in the large detuning case. If we assume e — > 







and Q being arbitrary positive real number, we have \I/ 0S ~ , which corresponds to the 

w 

carrier excitation. Both solutions above are also in good agreement with the solutions in a 
Fock state representation' 12 ! . However, from Eq.(8) we can know that, different from former 
solutions for the WER case, the probability amplitudes of up and down states of the ion are 
different and complicated in the case of the SER with large detunings. This is a new result 
which has not been obtained before by general approaches made in a Fock state or dressed 
state representation. 

For a more general comparison, we should first present the solutions under the RWA. 
As referred to in Sec. I, Eq.(l) can be treated to be a nonlinear JC model with the ap- 
proach proposed in Ref.[8] under the RWA. However, to make the comparison more easy 
and clear, we start our RWA treatment from Eq.(2). By performing a unitary transformation 
exp[-i(^o- 2 + T^r)t] on Eq.(2) with M — 1, 2, 3, • • •, we can obtain 

H M = (1- 2- M )a + a + g(a + a_ + aa + ) + g 2 (9) 

corresponding to the resonance conditions of = 2~~ M under the RWA, and the eigenenergy 
in a Fock state representation is of the form 

1 T] 2 1 



E± = (1 - 2"")(n + -) + -L ± -yV ( n + !) + (!- 2^) 2 (10) 

with n — 0, 1, 2, • • •. Meanwhile, performing a unitary transformation of exp[— i(^<r 2 + ^^)t] 
on Eq.(2) with K — 1, 2, 3, • • •, under the RWA resonance conditions of Q = K, we have 

H K = ^-^na z +g(a + a_ + aa + )+g 2 (11) 

whose eigenenergy in a Fock state representation is of the form 




(12) 



withn = 0, 1,2, 

As shown in the last section, although our method can in principle present all series 
solutions, we merely present the specific forms of the non-RWA solution for two simplest 
cases. Here we will use those two cases for a comparison with Eqs.(lO) and (12). For the 
case (1), we find that the two solutions for E = 1 ± e, along with two different restricted 
relations respectively, actually correspond to the same expression 



However, for the case (2), the situation is somewhat complicated. The direct algebra shown 
in Appendix presents us four solutions of E for this case, whereas the specific calculation 
shows that the solutions of Eqs.(A3) and (A4) are actually identical. 

In Figs. 2 and 3, the solutions under the RWA are compared with those for the non-RWA 
treatment in the case of f2 =0.5 and 3.0. We can find that, in the case of small values of f2 and 
rj, some of the solutions under the RWA can approach those without the RWA. It is physically 
reasonable because the RWA is only valid for the WER and LDL. However, although the 
solutions without the RWA are exact, only a few of these solutions can be obtained for a 
certain termination of the series solution. We can not obtain all the particular solutions of 
the system unless we make the termination of the series solution for almost infinite times up 
to the case with infinite series terms. In contrast, the RWA treatment can present general 
solutions, like Eqs.(10) and (12). Nevertheless, as it merely retains some resonance terms 
in the Hamiltonian, a specific solution under the RWA only corresponds to a specific value 
of f2. Moreover, the figures tell us that, only when the value of rj is not taken to be larger 
than 0.1, can we consider that the solutions under the RWA and without the RWA are in 
good agreement in the case of the WER. The RWA description is obviously inaccurate for 
the case beyond the LDL although is less than 1.0 in this case. When the value of f2 is 
larger, the difference between the RWA case and non-RWA one becomes great, as shown in 



E — - 




(13) 
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Figs. 3. It is difficult to find any correspondence for the two cases. However, there are some 
crossing points between the RWA solutions and non-RWA ones. It means at certain values 
of r), the solutions for eigenenergies in these two cases can be the same. In physics, it can 
be considered that the non-RWA solutions are some isolated ones, merely corresponding to 
those of up internal levels of the ion. One may not expect to obtain the general eigenenergies 
with our non-RWA treatment. However, the demonstration of above series solutions means 
the existence of the finite dimensional invariant subspace of some operators. Particularly 
from the simple forms of the eigenenergies, one may expect that the quantum KAM theorem 
should be applicable for the SER problem' 20 ! . 

Before ending our discussion of the usefulness of non-RWA solutions, it is also of interest 
to present a practical example of preparing the nonclassical motional states of the ion with 
our solutions. As referred to in Sec. I, there is a considerable interest in investigating the 
SER, particularly for rapidly preparing nonclassical motional states of the ion. In fact, 
with the results in the present paper, one can also generate some particular nonclassical 
states through suitable adjustment of the Rabi frequencies and proper measurements. For 
example, consider the simplest case, i.e., the case (1) in Sec II. For the case of the resonance 
(e = 0), we have E 1 = E' x with the same f2 and double degenerate eigenfunctions ^f(a) 
and ty'(a). Transforming our results into the Schrodinger representation will yield the 
states = exp (-iEtf + iuJia z t/2)UV^(a) and $'(*) = exp {-iE x t + iuj l a z t/2)UV^' (a), 
respectively. In the case of the LDL(g <C 1), the Rabi frequency f2 would be 2y/l — 4g 2 ~ 2. 
So we have c ~ ^zi&o, h ~ p^&o, c' Q ~ ^-T b 'o and b 'i ~ P^T fo o- B y controlling the 
evolution time t = Air/cui, and measuring the excited state of the ion, we can obtain a 
displaced even coherent state, a kind of Schrodinger-cat states 

$ M - ^ D (P)(\il >+\-^l>) = ^(N > +|0 >) (14) 

where the relation between the coherent state and Fock state has been used. As the Rabi 
frequency approaches 2, such a preparation is obviously outside the WER. 
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IV. CONCLUSION 



As e ga a n corresponds to the displaced Fock states \n,g >, the eigenf unctions of the 
system we obtained are infinite superposition of displaced Fock states. Therefore, it is 
understandable that we can not obtain similar analytical results by general Fock state ex- 
pansion. The obvious advantage of our approach is the possibility to obtain some relations 
between parameters of the system corresponding to certain nonclassical states by truncating 
the series expansion step by step. So as long as we can reach these conditions experimentally, 
the different nonclassical states predicted by our theory would be obtained. 

In summary, the Raman interaction between a trapped ultracold ion and two travelling 
wave lasers has been treated analytically in an interaction representation, without considera- 
tion of the RWA and the limitation for both the LDL and WER. Although the coherent state 
expansion technique was used previously for treating a similar problem, and the solutions 
we presented here are some isolated ones under certain conditions related to Q, rj and A, the 
present investigation is more general and exact so that we can study the ion trap problem 
in a wide range of parameters and compare our results with other approximate works or 
numerical solutions in this respect as we did in the present paper and for the cavity-atom 
problem ^ . More importantly, the Raman process has been widely used for the laser-cooling 
in current ion trap experiments, and the work outside the WER is attracting much interest 
of experimentalists. One possible way' 10 ' to experimentally realizing the SER is that the 
ion is first cooled within the LDL and under the WER, and then the trap frequency is 
decreased by opening the trap adiabatically so that the ratio of the Rabi frequency to the 
trap frequency is increased to a large number. Therefore, we believe that the difficulties in 
realizing the SER will soon be overcome, and our work would be helpful for any possible 
future exploration of the quantum properties of the ion-trap system outside the WER. 

The author sincerely thanks M.A.Kornberg for his carefully reading the manuscript. The 
work is partly supported by National Natural Science Foundation of China. 
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V. APPENDIX 



For obtaining a more specific expression of the eigenenergy for the case (2) in Sec. II, we 
can suppose X = e — g 2 for the case of z — g and Y = e + g 2 for z = —g. Direct algebra on 
Eq.(7) can present us following two second-order differential equations 

AX 2 + BX + C = (Al) 



and 



AY 2 - BY + C = (A2) 



with A = 8(1 - g 2 ), B = (3 + § )(2 + ft)(2 - § ) - 3(f - 2 + f + 2g 2 ) - 28g 2 - (3 + ft) (2 + 
§ - f - 2^ 2 ) and C = -[(f - 2 + f + 2^ 2 )(1 + Q)(2 - 0) + 20^ 2 - Qg 2 (^ + 2^ 2 )]. So we 
can obtain the eigenenergies 



^f = 2 + ^ + *— A (A3) 



for the case of z — g and 



1± o , , 5 ± VS 2 - 4AC 



^2=2 + ^ ^ (A4) 



for the case of z — — g. 
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Captions of the figures 



Fig.l Level scheme of the internal structure of the trapped ultracold ion, where \g ><-> 
|e > is dipole forbidden. 

Fig. 2 Variation of E with respect to rj in the case of Q = 0.5, where solid curves from 
the bottom to top are RWA solutions of Eq.(10), corresponding to n =0, 1, 2, 3, 4, 5 and 6, 
respectively, dot-dashed curve is the solution of Eq(13) and dashed curves are from Eq.(A3) 
or Eq.(A4) (because the solutions of Eqs.(A3) and (A4) are identical), (a) demonstrates 
E + , and (b) is for E~ . 

Fig. 3 Variation of E with respect to rj in the case of f2 = 3.0, where solid curves from 
the bottom to top are RWA solutions of Eq.(12) for the case of E + , corresponding to n =0, 
1, 2, 3, 4, 5 and 6, respectively, the dot-dashed curve is the solution of Eq(13) and dashed 
curves are from Eq.(A3) for cases of E ± respectively. 
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